typeset using JPSJ.sty <ver.l.Ob> 



One-Dimensional S 



1 Spin-Orbital Model with Uniaxial Single-Ion 
Anisotropy 



O 

o 
o 

Q 

(N 



X3 
O 

o 



> 
o 

O 



O 

o 



X 



Satoshi MiYASHlTA * and Norio Kawakami 

Department of Applied Physics, Osaka University, Suita, Osaka 565-0871, Japan 
(Received February 2, 2008) 

We investigate ground-state properties of a one-dimensional 5* = 1 spin-orbital model with or 
without uniaxial single-ion anisotropy. By means of the density matrix renormalization group 
method, we compute the ground-state energy, the magnetization curves and the correlation 
functions. We discuss how the ground-state properties depend on the two exchange couplings 
for orbital and spin sectors. The phase diagram obtained is compared with that for the S = 1/2 
model. We also address the effect of uniaxial single-ion anisotropy. 
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§1. Introduction 

The interplay of spin and orbital degrees of freedom 
has provided a variety of interesting phenomena in cor- 
related electron systems. To understand the role played 
by spin and orbital fluctuations, several different spin- 
orbital models have been studied extensively. As the 
simplest example among others, the one-dimensional 
(ID) spin-orbital model with SU(4) symmetry, which 
describes a spin-1/2 system with two-fold orbital de- 
generacy, has been investigated.^"*^ A slightly extended 
model, the ID spin-orbital model with SU(2)(8)SU(2) 
symmetry,'"^"*) has also been studied and the ground- 
state phase diagram has been established, which consists 
of a variety of phases including gapful/gapless spin and 
orbital phases, etc. 

Some transition metal oxides such as manganetes and 
vanadates, for which the Hund coupling plays an im- 
portant role, have higher spins with degenerate orbitals. 
This naturally motivates us to extend the above spin- 
orbital model to a model possessing higher spins. A spe- 
cific generalization of the model to the S = \ case has 
been proposed for vanadates, such as YVO3, and inves- 
tigated in detail. ^"^'^^ In particular, the competition of 
the orbital-valence-bond (OVB) solid phase and spin- 
ferromagnetic phase has been clarified, in accordance 
with some experimental findings in the low-temperature 
quasi-lD phase of YVO3. The realization of the OVB 
phase has been suggested in neutron diffraction experi- 
ments."'^"'-' 

Motivated by the above hot topics, we study an 
SU(2)(g)SU(2) extension of the S' = 1 spin-orbital model 
in ID. In contrast to the model proposed for YVOa,^-* 
our aim is to capture generic features inherent in the 
S = 1 spin-orbital model, which is to be compared with 
the S* = 1/2 model. ^'•^"^^ We also take into account the 
effects of single-ion anisotropy, which may play an im- 
portant role for S = 1 systems. We exploit the density 
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matrix renormalization group (DMRG) methodic) and 
investigate quantum phase transitions. By calculating 
the ground state energy and the spin/orbital correlation 
functions for a given strength of single-ion anisotropy, we 
obtain the phase diagram in the plane of two exchange- 
coupling constants for spin and orbital sectors. 

This paper is organized as follows. After a brief ex- 
planation of the model in the next section, we present 
the DMRG results in § 3 for the ground-state energy, 
the spin/orbital magnetization curves and the correla- 
tion functions, from which we determine the phase dia- 
gram. In § 4, we address the effects of uniaxial single-ion 
anisotropy. A brief summary is given in § 5. 

§2. Model 

We consider an = 1 extension of the ID spin-orbital 
model with uniaxial single- ion anisotropy D, which is 
characterized by two coupling constants A and B re- 
lated to the spin and orbital degrees of freedom. The 
Hamiltonian reads 



(S, • S,+i + A) T, • T,+i + - 



(1) 



where Si is an S" = 1 spin operator at the «-th site and T.^ 
is a T = 1/2 pseudo-spin operator acting on the doubly- 
degenerate orbital degrees of freedom. J controls the 
magnitude of the exchange couplings, which will be taken 
as the energy unit in the following discussions. The D- 
term represents uniaxial single-ion anisotropy, which has 
been discussed in detail so far in the Haldane spin chain 
systems. 

At a special point {A = B = 1 and D = 0), the sym- 
metry is enhanced to SU(2)(8)SU(2). It is known that the 
ground state in this case is the orbital liquid with a small 
spin gap, which is called the OVB solid state. ^"'^'^-' Ac- 
tually, the special model studied by KhaliuUin et al. for 
cubic vanadates^' '^^""''^-' includes the above special case, 
at which our model coincides with theirs. 
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The model (1) is regarded as a natural extension of 

the ID 5 = 1/2 spin-orbital model with SU(2)(8)SU(2) 
symmetry,^' for which the D-torm is absent, and A — > 
A/ A for the spin {S — 1/2) part. The comparison of 
these two models may allow us to clarify the role of spin 
fluctuations on the spin-orbital model. 

§3. Ground State Properties without Anisotropy 

In this section, we investigate the ground-state prop- 
erties of the Hamiltonian (1) with D = 0, and deter- 
mine the zero-temperature phase diagram. Wc first no- 
tice that the spin- (orbital-) ferromagnetic state should 
be the ground state for B — > — oo {A —> — oo). Therefore, 
when wc calculate the groimd-state energy of the spin- 
(orbital-)ferromagnetic state, we can fix (S, • Sj+i) = 1 
((T, • T,+i) = 1/4). 



the exact ground-state energy for the isotropic T = 

Heisenberg chain. '^^^ On the other hand, for B - 
the I-IE boundary approaches A = \eg^^\, where e 



= 1/2 

^ oo. 



—1.401484 is the ground-state energy for the S = 1 

Haldane chain. This asymptotic analysis is consistent 
with our phase diagram obtained here numerically. 
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Fig. 2. Phase diagram of the ID 5 = 1 spin-orbital model in the 
A-B plane without uniaxial single-ion anisotropy. The phase I 
is the OVB phase (see text). In the phase 1, the spin part is in 
the fully polarized state while the orbital part forms the gaplcss 
antifcrromagnctic T = 1/2 Heisenberg chain. On the other hand, 
in the phase IH, the orbital sector is in the ferromagnetic phase 
while the spin sector is in the S = 1 Haldane phase. Both of the 
spin and orbital parts are in the fully polarized ferromagnetic 
states in the phase IV. All the transitions are of first order. 



Fig. 1. The energy per site as a function of the exchange param- 
eter B for A = D = 0; the ground-state energy obtained by 
the DMRG is shown by triangles. Wc also plot the energy com- 
puted for the spin-fcrromagnctic state (circles) and the orbital- 
ferromagnetic state (squares). We clearly sec two cusp struc- 
tures, indicating first-order phase transitions. The intermediate 
region is identified as the OVB solid state. 



In Fig. 1 , the energy obtained by the DMRG is shown 
as a function of B with keeping A = fixed. We 
find two first-order transition points; Bc\ — —0.09 and 
Bc2 — 0.41. As clearly seen in Fig. 1, the system favors 
the spin-fcrromagnctic state for small B (< Bci), but 
stabilizes the OVB solid state'^' for the intermediate re- 
gion {Bel < B < Bc2)- For large B (> Sc2), the ground 
state is in the orbital-ferromagnetic phase. 

By repeating similar estimations of the critical points 
for other choices of A, we determine the ground-state 
phase diagram for D = in the A-B plane, which is 
shown in Fig. 2. The phase diagram consists of four 
phases; the OVB phase (I), the spin- ferromagnetic phase 
(I), the orbital- ferromagnetic phase (IH), and the spin- 
ferro/orbital- ferromagnetic phase (IV). We can estimate 
the phase boundary exactly in several limiting cases. 
First, by dropping the spin (orbital) degrees of freedom 
when the spin- (orbital-) ferromagnetic ground state is 
stabilized, we immediately find that the phase bound- 
aries of E-IV and IE- IV are exactly given hy A = —1 and 
B = —1, respectively. We can also determine the ex- 
act asymptotic behavior of the boundary. The boundary 
between the phases I and I approaches B =- 4|eg 
for A — > oo, where Sg^"^^ 



■log2 + 1/4 
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Fig. 3. Spin (upper panel) and orbital (lower panel) correlation 
functions in the phase I along the A = B{= x) line in the OVB 
phase. The site index i is even (odd) for open circles (solid 
triangles). The even-odd i-dependence reflects a dimer property 
of the OVB phase. With increasing x, the dimerization gets 
strong and the size of a gap becomes large. 



Although characteristic properties of the phases E, 
in and IV can be directly deduced from those of the 
S = 1/2, 1 chains, there appear nontrivial properties due 
to the interplay of spins and orbitals in the OVB phase I. 
We thus calculate the short-range spin and orbital cor- 
relation functions by taking the line A = B[= x) as 
representative parameters in the OVB phase. 

The results are shown in Fig. 3. As x increases, the 
even-odd i-dependence of the correlation functions gets 
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Fig. 4. Spin (a) and orbital (b) magnetization curves at A = B = 
1 (solid line) and 0.5 (dashed line). We calculate the spin (or- 
bital) magnetization by setting ho = 0{hs = 0). Note that both 
of the spin and orbital magnetizations exhibit typical behavior of 
gapful systems, although the orbital gap is much larger than the 
spin gap. For higher fields, the spin and orbital magnetizations 
show the first-order phase transition to the fully polarized state, 
which is signaled by a discontinuous jump. 



stronger both for spin and orbital sectors, implying that 
the system favors the dimerization. At a; = 1, the orbital 
sector forms nearly perfect dimer singlets, i.e. (T^ • 
Tj+i) ~ —3/4 for even i, while it is almost zero for odd 
i. On the other hand, two adjacent spins are almost 
parallel for even i, i.e. (Sj • S^+i) ^ 1, while for odd 
i, there are weak antiferromagnetic correlations. These 
properties clearly characterize the dimer-like nature of 
the OVB phase.^' It is remarkable that even if x 
decreases from unity, the orbital correlation (Tj • Tj+i) 
for odd i stays very small, indicating that the orbital 
sector can be regarded as an assembly of approximately 
independent dimers irrespective of the values of x. 

In order to observe the properties of the phase I in 
more detail, wc compute the spin/orbital magnetization. 
To this end, we add the following terms to the Hamilto- 
nian (1), 

Heo.^-h,J2^^-''oJ2T;■, (2) 
i i 

where and ho are external fields conjugate to the 
spin and orbital magnetizations. The calculated mag- 
netization curves are shown in Fig. 4 for two typical 
choices of the exchange couplings. Let us first observe 
the spin magnetization for x = A = B = 1 shown in 
Fig. 4(a). The magnetization (Sf) is zero in small fields, 
in accordance with the existence of a spin gap in the 
OVB phase. After the spin gap disappears at very low 
fields, the magnetization increases smoothly, and jumps 
at the second critical point (e.g. hg ~ 0.122 for x = 1), 
driving the system to the spin-ferromagnetic state via a 
first-order transition. Similarly, the orbital magnetiza- 
tion shown in Fig. 4(b) indicates that the OVB state 
with {Tf) = is favored in small fields, and then the 



orbital magnetization (Tf) develops gradually beyond 
a critical field corresponding to the orbital gap. The 
system further undergoes a first-order transition to the 
orbital-ferromagnetic state. Comparing the spin and or- 
bital gaps, we notice that the orbital gap is much larger 
than the spin gap,^^'^'^'' consistent with the results for 
correlation functions, where the orbital-dimerization is 
quite strong. In contrast to the low-field behavior, both 
spin and orbital degrees of freedom change their charac- 
ters at the first-order transition points. Note that the 
resulting high-field phase with the fully spin (orbital) 
polarized state is the same as the phase II (IE). 
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Fig. 5. Phase diagram of the S = 1/2 spin-orbital model."''"**' 
Phase 1 is a gapful phase with doubly degenerate ground states 
which form alternating spin and orbital singlets. In the phase I, 
the spin degrees of freedom are in the fully polarized ferromag- 
netic state while the orbitals cire in the antiferromagnetic ground 
state and vice versa in the phase M. Both spin and orbital de- 
grees of freedom are in the fully polarized ferromagnetic states 
in the phase IV. The phase V is a gapless phase including an 
integrable SU(4) point (A = B = l). 



Before closing this section, we compare the present 
phase diagram with that for the 5 = 1/2 spin-orbital 
model, ^"^^ which is sketched in Fig. 5. The corre- 
sponding spin-orbital Hamiltonian is given by (1) by 
putting S = 1/2 and replacing A with A /A for the 
spin sector. Overall features are similar both in 5 = 
1 and 1/2 models: the phases I, H and IV are re- 
spectively characterized by the spin-ferro/orbital-singlet 
state, the spin-singlet/orbital- ferro state and the spin- 
ferro/orbital-ferro state. The phase I also exhibits sim- 
ilar dimer-like properties in both cases. A remarkable 
difference is the gapless phase V realized in the 5=1/2 
model around the SU(4) point. This phase, which is 
stabilized by enhanced quantum fluctuations both in 
spin and orbital degrees of freedom, disappears in the 
5=1 model. We think that the phase V is inherent 
in the 5=1/2 case, and a higher-spin extension of the 
SU(2)(8)SU(2) model may have a phase diagram similar 
to the present 5=1 case. 

In contrast to the 5 = 1/2 model, the single ion 
anisotropy plays an important role for the 5=1 model, 
which may provide a rich phase diagram. This problem 
is addressed in the following. 
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§4. Effects of Uniaxial Single- Ion Anisotropy 

4.1 D<0 

Let us now consider the effects of uniaxial single-ion 
anisotropy. We start with the ground-state properties 
for the case of _D < 0. 

Recall that in the phases I and IV in Fig. 2, the spin 
sector is in the fully-polarized state, so that the nature of 
these phases may not be altered upon the introduction of 
D <0, except that the direction of spin ordering is now 
fixed to the 2;-axis. On the other hand, the negative D 
suppresses quantum fluctuations in the spin-gap phases I 
and n, resulting in the antiferromagnetic spin order with 
Z2 symmetry breaking. We refer to the corresponding 
new ordered phases as the phase and IHm, respectively. 

Note that different types of spin order should emerge 
in the phases and Mm ■ The spin-antiferromagnetic or- 
der has the 4-sitc period (up-up-down-down alignment) 
in the phase I„i, and the 2-site period (up-down- up-down 
alignment) in the phase Wm, as schematically shown in 
Fig. 6. This difference comes from the ground-state 
properties without anisotropy. Namely, in the OVB 
phase I, spins favor a ferromagnetic configuration in a 
orbital-dimer singlet and an antiferromagnetic configu- 
ration between adjacent orbital singlets, so that the spin 
sector is described by an effective ferro-antiferromagnetic 
bond-alternating c;liain. On the other hand, in the phase 
Mm, the spin sector forms the S = I Haldane state at 
D = 0, so that it naturally leads to an ordinary antifer- 
romagnetic order upon the introduction of D < 0. 



• phase I (i=odd) 
o phase I (i=even) 
□ phase III 
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Fig. 6. Schematic description of spin orders for the phases Im and 
M^. The label S (T) denotes the spin (orbital) part, respectively. 



The above tendency to spin ordering is indeed ob- 
served in the nearest-neighbor correlation functions 
shown in Fig. 7 as a function of D{< 0). We exploit 
two typical values of A and which correspond to the 
parameters for the phases I and IH at = 0. Even in 
the presence of D, both spin and orbital correlation func- 
tions show (do not show) the even-odd i-dependence for 
the parameters corresponding to the phase 1^ (JS.m)- It 
should be also noticed that the orbital correlations are 
almost independent of _D, as should be naively expected: 
almost isolated orbital-dimer states are favored in the 
phase Im, while the orbital sector is always in the fully 
polarized state, (T^ ■ T^+i) = 1/4, in the phase Wm ■ 

As mentioned above, the spin-antiferromagnetic order 
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Fig. 7. Spin (upper and middle panels) and orbital (lower panel) 
correlation functions as a function of D; A = B = 1 (solid circles 
: i = odd and open circles ; i = even) and A = —1, B = 1 (open 
squares). 



with the 4 (2)-site period should be stabilized beyond a 
certain critical value of jZ?! in the phase Im (Hm). There- 
fore, spin fluctuations are suppressed and thus the cor- 
relation (5? • ^ ±1 while {S+ ■ S-^i) ^ as \D\ 
increases. We now determine the phase transition point 
for IH— >lEm. Note that the orbital sector always orders 
ferromagnetically, so that the Hamiltonian (1) for the 
spin sector is reduced to 

noF = J'Y^{Si-Si+i + A) + DY^{Stf, (3) 

i i 

with J' = J(l -|- B)/4. In this case, the competition of 
the single-ion anisotropy D and the effective exchangc;- 
coupling J' results in the following transition, which 
is independent of A. According to the results for the 
isotropic S = 1 chain, the ground state is in the 

Haldane phase for —0.35 J' <D< J', and in the Neel phase 
for D < — 0.35J'. Therefore, the critical value is given 
by Dc = —0.175 J in Fig. 7, which separates the Haldane 
phase and Neel phase for the spin sector. 

It is more difficult to determine the boundary between 
I and Im only from our data for the correlation functions 
in Fig. 7, partially because the spin gap is very small 
in the phase I. Nevertheless, by taking into account that 
the spin gap is much smaller than the orbital gap (Fig. 
4), we reasonably expect that the corresponding critical 
value of D, which separates I and in Fig. 7, may 
be quite small, say, less than —0.1. The precise critical 
value should be obtained by exploiting improved numeri- 
cal methods, which we wish to leave for the future study. 

We show the phase diagram in the A-B plane for 
£> < in Fig. 8. Since the transition from the OVB 
phase to the other phases is of first order, we deter- 
mine the corK^sponding boundary by comparing three- 
types of the energy among the OVB solid state, the spin- 
ferromagnetic state and the orbital-ferromagnetic state. 
As mentioned above, the phase boundary between the 
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phases M and Mr, 



is given by 
Br ^ -HAD 



1, 



(4) 



which is independent of A. For example, Be — 0.14 
for D = —0.1, as shown in Fig. 8(a). In the case of 
D = —1.0 in Fig. 8(b), the phase Mm dominates the 
phase M in the displayed region, because the boundary 
between these two phases is at Be — 10.4 for D = —1.0. 
Note that the orbital sector orders ferromagnetically for 
all D in the phases n. Mm and IV. As mentioned above, 
the spin state in the OVB phase I is sensitive to D, so 
that the ordered phase Im dominates the phase I both 
for D = —0.1 and —1.0 in the region shown in Fig. 8 
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Fig. 8. Phase diagrams for D < 0: (a)D = —0.1 and (b)D = 
— 1.0. The phases labeled as I, II and IV arc the same as in Fig. 
2. Im and III,„ are spin-ordered phases (see text) . The solid lines 
except for the M-Mm boundary represent first-order transitions. 
For reference, the phase boundary at D = is shown by the 
dashed lines. 



4.2 D>Q 

Let us now move to the positive D case. We start 
with the phases IH and IV, for which the phase transition 
points are easily determined, since the orbital sector is 
in the fully polarized state. In the phase M, the so-called 
large- Z) phase is stabilized in the spin sector for D > 
J' = J(l + B)/4; i.e. for -1 < B < B^ {B^ = -1 + 
4:0), the large-D phase appears, while the Haldane phase 
persists for B > B^- On the other hand, in the phase IV, 
the XY phase emerges^^'^^'^^-* in the region —1 — AD < 
i? < — 1 in the presence of the D-term. 

It turns out that the boundary between the phases I 
and II is quite sensitive to positive D, and exhibits some- 
what complicated features. This is contrasted to the 
robust first-order transition between the phases I and IE, 
which hardly changes its character as far as D is small. 
To make the above point clear, we observe the behavior 
of the correlation functions by choosing D = 0.01 with 
two typical values of ^4 = and 0.5. In Fig. 9 (a), the 
S-dependence of the spin and orbital correlation func- 
tions are shown for A = 0. For small B, the ground 
state is in the phase I, so that both spin and orbital 
correlation functions are spatially uniform without the 



even-odd i-dependence. With increasing B, a continu- 
ous quantum phase transition occurs from the phase I 
to I at i3c ~ —0.1, in contrast to the first-order transi- 
tion at £> = 0, as seen in Fig. 9(a). Beyond the critical 

value Be, the spin and orbital correlation functions ex- 
hibit the even-odd i-dependence, reflecting dimer proper- 
ties of the OVB solid state. Unfortunately, we cannot fig- 
ure out whether this continuous transition is Berezinskii- 
Kosterlitz-Thouless (BKT) type^^"^'"') or not only from 
the results shown in Fig. 9, because an extremely large 
system size is necessary to obtain the sensible results 
in our spin-orbital model with a tiny spin gap. Never- 
theless, we can discuss some characteristic properties of 
each phase from the correlation functions. First, we note 
that the spin sector in the phase II may be changed to 
the XY phase upon introduction of D.18,21,22) ^jjjg ^.^^^ 
be inferred from the spin correlation functions shown in 
Fig. 9, which indeed show the preference of the XY 
phase: (S+S-^^) ~ 1 and (5'fS'f+i) - 0. This obser- 
vation is also supported by the following consideration. 
Note first that the orbital correlation function is not af- 
fected by D in the phase H, and takes the constant value 
(T, • T,+i) - e^^/^' characteristic of the T = 1/2 orbital 
chain, so that we can focus on the remaining S — 1 spin 
part. According to the isotropic S = 1 ferromagnetic 
chain,^*'^^'^^' the first-order phase transition from the 
spin-ferromagnetic phase to the XY phase occurs upon 
introducing D, which is consistent with our numerical 
results for the correlation functions. 



(a) A=0.0 

1(^0000000000^, 



• i=odd 
o i=even 



(b) A=0.5 



+ 

CO 





-1 

0.4 
0.2 


-0.2 
-0.4 

02 

-0.2 



°oooooi 



••••••<l 



liooooooooaBtooooooooo; 



-I - 

0.4^ 
0.2- 



- -0.2- 

- -0.4- 



' ^ 0.2 





" 78088880880(800000000^ "°45O00O0OO0(D 

0.6k - -0.6 - n 



-0.1 
B 



0.15 



(gPOOOOOOOO^ 

0.2 0.25 0.3 

B 



Fig. 9. Spin (upper and middle panels) and orbital (lower panel) 
correlation functions for D = 0.01 as a function of B ((a) 
A = 0.0, (b) A = 0.5). The site index i is even (open circle) 
or odd (solid circle). The even-odd j-dependence reflects dimer 
properties of the phase I. 



The above continuous transition is in contrast to that 
for A = 0.5 shown in Fig. 9 (b). where the system ex- 
hibits the first-order phase transition accompanied by 
a clear jump at Be ~ 0.2 both in the spin and orbital 
correlation functions. This characteristic behavior is es- 
sentially the same as that for Z? = 0. 

By examining the correlation functions for other 
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choices of A, we end up with a phase diagram expected 
for D = 0.01, which is shown in Fig. 10. The phase 
boundary between the phases I and H is cither continu- 
ous or first-order depending on the value of A: the tran- 
sition is continuous (first-order) for A < 0.1 {A > 0.1). 
We have checked that the critical value of A, which sep- 
arates these two types of transitions, increases monoton- 
ically with the increase of D. For example, the critical 
value is A = 0.5 for D = 0.012. For smaU D, the first- 
order transition dominates most part of the I-I phase 
boundary. 

Here, an important question arises: what really causes 
the change in the nature of the transitions, first-order 
or continuous, on the I- II phase boundary. To address 
this question, we should go into the detail of the phase 
I. Since the orbital part always forms the dimerized 
state in the phase I even for finite D, we again focus 
on the spin sector. Recall here that the spin sector in 
the OVB phase I forms the state similar to the 5* = 2 
Haldane-gap state (same universality class). For- 
tunately, the ID 5* = 2 Haldane system with uniax- 
ial single- ion anisotropy was already investigated.^^' ^'^^ 
Schollwock numerically checked the disappearance of the 
spin gap at a very small critical value of D (~ 0.02J), 
and the system enters the " intermediate- £) phase" , which 
was claimed to be equivalent with the XY-phase.^''^ Im- 
mediately after the above analysis, Oshikawa found that 
the behavior of the string correlation in z-axis devi- 
ates from a simple power-law, as a result, he suggested 
that this intermediate-!) phase is not truly gapless XY- 
phase, but somewhat modified one with some anomalous 
properties in the 2;-components.^^) Anyway, the message 
from the above analysis is that the introduction of small 
positive D changes the Haldane phase to the so-called 
intermediate-D phase. In our model, the size of the spin 
gap in the phase I gets small as A — > — 1 and B —1 
(Fig. 4). Therefore, we think that the spin sector in 
the OVB phase is changed to the " intermediate- £> type 
phase" in the region with small A and B, while the orig- 
inal spin-gap state in the OVB phase can still persist for 
larger A and B. This may cause the change in the na- 
ture of the phase transition on the I-I phase boundary. 
Unfortunately, our numerical calculation is not powerful 
enough to draw a definite conclusion about the existence 
of the "intermediate-!? type phase". We would like to 
clarify the detail in the future work. 

§5. Summary 

We have investigated quantum phase transitions of the 
ID 5 = 1 spin-orbital model with uniaxial single-ion 
anisotropy. By means of the DMRG method, we have 
calculated the ground-state energy, the spin/orbital mag- 
netizations, and the correlation functions, from which we 
have determined the zero-temperature phase diagram. 

In the absence of the anisotropy, there appear four 
phases including the OVB phase characteristic of the S = 
1 spin-orbital system.^) This phase has a small spin gap 
and a fairly large orbital gap. In comparison with the 
5=1/2 spin-orbital model, we have found that most 
of the phase diagram shows similar properties, but the 
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Fig. 10. Phase diagram for the model with D = 0.01. The solid 
lines (with dots) indicate the first-order transition. Concern- 
ing the phase boundary between I and I, there is a continuous 
transition (broken line) for A < 0.1, and a first-order transition 
(solid line with dots) for A > 0.1. The intermediate-D phase 
is expected to emerge in the spin sector of the OVB phase for 
small A and B. Note also that around the boundary oi B = — 1, 
there appear two new spin phases: the large- 1? phase (XY spin 
phase) in the tiny region -1 < B < -0.96 (-1.04 < B < -1), 
although not shown explicitly in the figure (see text). 



phase V inherent in the S* = 1/2 model disappears in the 
S = 1 case. We believe that a higher spin extension of 
the model should have the phase diagram similar to that 
of the present 5 = 1 model. 

The introduction of uniaxial single-ion anisotropy D 
gives rise to some interesting aspects. For negative D, 
two different types of magnetic orders appear, depend- 
ing on whether the system at £) = is in the orbital- 
ferromagnetic phase or the OVB phase. Since the spin 
gap for the OVB phase is much smaller than that for the 
orbital-ferromagnetic phase, the spin order in the OVB 
phase is induced even in the small \D\ region. 

On the other hand, it has turned out that the situation 
is much more subtle in the positive D case. In particu- 
lar, the phase transition between the OVB phase I and 
the spin-ferromagnetic phase H exhibits somewhat com- 
plicated feature: the nature of transition changes from 
first-order to continuous one, depending sensitively on 
the value of D. Although we have not been able to ob- 
tain the precise phase diagram, some characteristic prop- 
erties of each phase have been discussed on the basis 
of the correlation functions. These nontrivial properties 
originate from interplay of the spin and orbital degrees 
of freedom, which may exemplify interesting aspects of 
the spin-orbital systems. Further investigations should 
be done in the future study, especially for the positive D 
case. 
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